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Incorporating a newly improved isospin- and momentum-dependent interaction in the isospin- 
dependent Boltzmann-Uehling-Uhlenbeck transport model IBUU11, we have investigated relative 
effects of the density dependence of nuclear symmetry energy E sy m{p) and the neutron-proton effec¬ 
tive mass splitting — m y on the neutron/proton ratio of free nucleons and those in light clusters. 

It is found that the m£ — has a relatively stronger effect than the E aym {p ) and the assump¬ 
tion of mj) < m,p leads to a higher neutron/proton ratio. Moreover, this finding is independent of 
the in-medium nucleon-nucleon cross sections used. However, results of our calculations using the 
E ayrrl (p) and m* — m* both within their current uncertainty ranges are all too low compared to the 
recent NSCL/MSU double neutron/proton ratio data from central 124 Sn+ 124 Sn and 112 Sn+ 112 Sn 
collisions at 50 and 120 MeV/u, thus calling for new mechanisms to explain the puzzlingly high 
neutron/proton ratio observed in the experiments. 
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To pin down the density dependence of nuclear sym¬ 
metry energy E sym (p ) has long been a major challenge 
for both nuclear physics and astrophysics. While much 
progress has been made in the past decade, many in¬ 
teresting issues remain to be resolved 00 . A larger 
symmetry energy generally corresponds to a more re¬ 
pulsive (attractive) underlying symmetry potential U syrn 
for neutrons (protons) in neutron-rich nuclear matter as 
they are linearly proportional to each other according 
to the Hugenholtz-Van Hove theorem [5] or the Bruck¬ 
ner theory 10, see, e.g., Refs. (81177)1 for the explicit 
relationship between E sym (p) and U sym . On the other 
hand, the in-medium nucleon effective mass describes to 
the first order effects due to the non-locality of the un¬ 
derlying nuclear interactions and the Pauli exchange ef¬ 
fects in many-fermion systems jll|. It can be calculated 
from the momentum dependence of the single-particle 
potential in non-relativistic models or the Schrodinger- 
equivalent potential in relativistic models. The nucleon 
effective mass is related to many interesting problems in 
both nuclear physics and astrophysics [12l - ll5l ] . It has fur¬ 
ther been found that the neutrons and protons may have 
different effective masses in neutron-rich matter due to 
the momentum dependence of the symmetry (isovector) 
potential. However, calculations within different models 
using various interactions, e.g., the Brueckner-Hartree- 
Fock approach [Ij| , the relativistic mean-field model 0, 
and the Skyrme-Hartree-Fock calculation [l8|, [n|, pre¬ 
dict rather different values for the neutron-proton effec- 
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tive mass splitting m*_ p = m* — m*. Thus, currently 
there is no consensus as to whether the m*_ p is nega¬ 
tive, zero, or positive. However, the value of m*_ p af¬ 
fects significantly isospin-sensitive observables in heavy- 
ion collisions 1 2(1 - 1271 as well as thermal and transport 
properties of neutron-rich matter [28H30I. It also has im¬ 
portant ramifications in astrophysics 1311 . For instance, 
the equilibrium neutron/proton ratio in the primordial 
nucleosynthesis is determined by ( n/p) eq = in 

the early (> lms ) un iverse when the temperature T was 
high (> 3 MeV) [32} . 

Recently, analyses of the free neutron/proton double 
ratio from central 124 Sn+ 124 Sn and 112 Sn+ 112 Sn colli¬ 
sions at 50 and 120 MeV/u at the NSCL/MSU seem to 
indicate that protons have a slightly larger effective mass 
than neutrons based on comparisons with calculations 
within an improved quantum molecular dynamics model 
using Skyrme interactions (33j. The preferred Skyrme 
interaction SLy4 [34} is being widely used in describing 
the ground state properties and excitations of neutron- 
rich nuclei. However, the applicability of Skyrme inter¬ 
actions is restricted, in general, by the nuclear struc¬ 
ture calculations and the small amplitude nuclear mo¬ 
tions due to the incorrect energy behavior of nucleon- 
nucleus isovector optical potential in comparison with 
that extracted from experiments [l(| f77[T~iT | . This situa¬ 
tion clearly calls for more theoretical studies with differ¬ 
ent transport models and examinations of various model 
ingredients. In fact, it was known that the isospin- 
dependent Boltzmann-Uehling-Uhlenbeck (IBUU) trans¬ 
port model using a momentum-dependent symmetry po¬ 
tential corresponding to m* > m* [42} under-predicts 
the old NSCL double neutron/proton data [43[ no matter 
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how the density-dependence of symmetry energy and the 
in-medium nucleon-nucleon cross sections are adjusted. 
Thus, to understand the puzzlingly high neutron/proton 
double ratio, the relative effects of the symmetry energy 
and the neutron-proton effective mass splitting should 
be studied within the same model. We have recently 
incorporated a newly improved i sosp in- and momentum- 
dependent interaction (ImMDI) [29| based on the Gogny 
force in the IBUU11 transport model. The original 
Gogny interaction is known to give an asymptotic value 
of the isoscalar potential at saturation density less than 
that extracted from optical model analyses of nucleon- 
nucleus scattering data. We removed this drawback and 
introduced a new parameter to adjust easily the values 
of the m* _ p besides those varyingly the magnitude and 
density dependence of the E sym (p). In this Brief Report, 
we revisit the neutron/proton puzzle. We found that in¬ 
deed the m* — m* has a relatively stronger effect than 
the E sym (p ) and the assumption of m* < m* leads to 
a higher neutron/proton ratio. However, the puzzle re¬ 
mains as our calculations using the E sym (p) and m* — m* 
both within their current uncertainty ranges still under¬ 
predict significantly the NSCL/MSU data. 

The ImMDI interaction is developed from the MDI 
interaction, which has a similar functional form as the 
Gogny effective interaction while replacing the Gaussian- 
type finite-range term with a Yukawa form [44j, [45| . The 
potential energy density in asymmetric nuclear matter 
from the MDI interaction or the ImMDI interaction is 
expressed as [3 
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In the above, p n and p p are respectively the neutron and 
proton density, and p = p n + p p is the total density, po is 
the saturation density, and S = (p n — p p )/p is the isospin 
asymmetry, r = 1(—1) denoting neutrons (protons) is 
the isospin index, and f T {r,p) is the phase-space distri¬ 
bution function. The single-particle potential from the 
mean-field approximation depends on the density p and 
isospin asymmetry 5 of the nuclear medium as well as the 
isospin t and momentum p of the nucleon [44| 
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Comparing with the MDI interaction, the ImMDI in¬ 
teraction has been improved mainly in two aspects [29;]. 
First, the high-momentum part of the nucleon isoscalar 
mean-field potential has been refitted to reproduce the 
optical potential extracted from the proton-nucleus scat¬ 
tering experimental data up to the nucleon kinetic energy 
of about 1 GeV. Second, besides the x parameter, which 
was previously used to mimic the density dependence of 
the symmetry energy by adjusting the relative contri¬ 
butions of different spin-isospin channels of the density- 
dependent interaction, another two parameters y and z 
are introduced to vary respectively the isospin splitting 
of the nucleon effective mass and the value of the symme¬ 
try energy at saturation density. The parameters y and 
z enter the functional form through 
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where p/o = ^(37r 2 po/2) 1 ^ 3 is the nucleon Fermi momen¬ 
tum in symmetric matter at saturation density. The val¬ 
ues of Aq, C u o, Cio, B, cr, and A are fixed by six empirical 
constraints at x = 0, y = 0, and z = 0, i.e., the saturation 
density po = 0.16 fu+ 3 , the binding energy Eo{po) = 
— 16 MeV, the incompressibility Kq = 230 MeV, the 
isoscalar effective mass m* 0 = 0.7 m, the symmetry en¬ 
ergy E sym (po) = 32.5 MeV, and the isoscalar potential at 
infinitely large momentum 17o,oo = 75 MeV, and the val¬ 
ues of the corresponding parameters are Ao = —66.963 
MeV, C u0 = -99.7017 MeV, C i0 = -60.4860 MeV, 
B = 141.963 MeV, a = 1.26521, and A = 2.42401p /o . 

The ImMDI interaction provides us with more flexi¬ 
bility to investigate the detailed isovector properties of 
nuclear interaction. In the present work, we set z = 0 
and vary the values of x and y to study effects of the 
symmetry energy and the neutron-proton effective mass 
splitting. The nucleon effective mass is defined as 


m* A m dU T \ 1 

m \ p dp J 


(7) 


and it generally depends on the density and isospin asym¬ 
metry of the medium as well as the isospin and momen¬ 
tum of the nucleon. The density dependence of the sym¬ 
metry energy E sym and the relative neutron-proton effec¬ 
tive mass splitting are displayed in Fig. [l] with different 
values of x and y. As discussed and shown in Ref. [29| . 
x affects only E sym while y affects both E sym and the 
m*_ p . In the following, we select several special sets of 
parameters to examine the relative effects of E sym and 
m*_ p . With (x = 0, y = —115 MeV) and (x = 1, y = 115 
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MeV), the symmetry energy is almost the same while the 
relative neutron-proton effective mass splittings are op¬ 
posite in sign. On the other hand, with (x = 0, y = —115 
MeV) and (x = 1, y = —115 MeV) the relative neutron- 
proton effective mass splitting is the same while the latter 
gives a softer symmetry energy. We can thus study the 
effect of the isospin splitting of nucleon effective mass by 
comparing the results from the former two parameters 
sets while investigate that of the symmetry energy by 
comparing the results from the latter two sets. We note 
that the current uncertainty range of the slope parame¬ 
ter L of E sym is about 50 ± 20 MeV [46[ , which is not 
quite different from the chosen range (10, 60) MeV in the 
present study. 




p/p 0 p (MeV/c) 


FIG. 1: (Color online) Density dependence of the symme¬ 
try energy (a) and momentum dependence of the relative 
neutron-proton effective mass splitting in nuclear matter of 
the density p = po and isospin asymmetry 8 = 0.2 (b) from 
the ImMDI interaction with different values of x and y. 


The ImMDI interaction with parameter sets described 
above was implemented in the IBUU11 model Q. In our 
calculation, 200 test particles per nucleon are used and 
about 20,000 events are generated for each beam energy 
and impact parameter. The initial density distribution is 
generated from the Skyrme-Hartree-Fock calculation us¬ 
ing the MSL0 force [471 > and the initial nucleon momenta 
are generated from local Thomas-Fermi approximation. 

The neutron/proton ratio in collisions induced by 
neutron-rich nuclei was firstly used as a probe of the 
symmetry energy [48j]. Later, the double neutron/proton 
ratio of nucleon emission for two collision systems with 
isotopes of different total isospin asymmetries was in¬ 
troduced to reduce systematically the influence of the 
Coulomb force and the poor efficiency of detecting neu¬ 
trons [42j, [49[ . To explore the effects of both the sym¬ 
metry energy and the isospin splitting of nucleon ef¬ 
fective mass on double neutron/proton ratio within the 
ImMDI and IBUU framework, we generate events for 



FIG. 2: (Color online) The coalescence invariant double 
neutron/proton ratios DR(n/p) in 124 Sn+ 124 Sn collisions to 
112 Sn+ 112 Sn collisions as a function of nucleon center-of- 
mass energy at beam energies of 120 (a) and 50 MeV/u (b) 
with the impact parameter b = 2 frn and the angular gate 
70° < 9cm < 110°. The NSCL data are from Ref. l33ll . 


112 Sn+ 112 Sn and 124 Sn+ 124 Sn collisions at beam ener¬ 
gies of 120 and 50 MeV/u. Similar to the treatment in 
Ref. we stop the evolution at t = 150 fm/c when 
the interaction becomes negligible, and identify nucleons 
and clusters based on the final nucleon phase-space dis¬ 
tribution, i.e., two nucleons are within one cluster if their 
spatial distance is closer than A r = 3 fm and their mo¬ 
mentum distance is smaller than A p = 300 MeV/c. We 
notice that our final results are not sensitive to the vari¬ 
ation of these coalescence parameters within about 30% 
of the above values. 

The coalescence invariant yield is constructed by con¬ 
sidering both free nucleons and those bound in light clus¬ 
ters including deuterons, tritons, 3 He, and 4 He, and the 
angular gate is chosen to be 70° < 9 cm < 110°, as in 
the experimental analysis. The impact parameter is set 
to be b = 2 fm to mimic the centrality in the exper¬ 
iments. The double neutron/proton ratio DR(n/p) in 
124 Sn+ 124 Sn and 112 Sn+ 112 Sn collisions is defined as 


DR(n/p) 


\Y(n)/Y(p)\ 124 Sn+ 124 Sn 
[V(n)/V(p)]ii2 S n+ii2Sn 


( 8 ) 


Since the yield neutron/proton ratio Y(n)/Y(p) in 
124 Sn+ 124 Sn collisions is larger than that in 112 Sn+ 112 Sn 
collisions, the DR(n/p) is always larger than 1. Figured 
shows the DR(n/p) as a function of nucleon center-of- 
mass energy at beam energies of 120 and 50 MeV/u. At 
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50 MeV/u, it is seen that the DR(n/p) for high-energy 
nucleons is slightly larger for a softer symmetry energy 
(L = 10 MeV) consistent with the finding in Ref. [H}. It 
is noteworthy that a stiffer symmetry energy can lead to 
a larger DR(n/p) at beam energies as high as 400 MeV/u, 
when the symmetry energy at suprasaturation densities 
becomes important. At the beam energy of 120 MeV/u 
shown in panel (a) of Fig. [2] the DR(n/p) is rather in¬ 
sensitive to the stiffness of the symmetry energy. On 
the other hand, the neutron-proton effective mass split¬ 
ting has a more appreciable effect on the DR(n/p). It is 
seen that a negative m*_ p results in a larger DR(n/p) at 
higher nucleon energies, consistent with our expectation. 
However, even for the two extreme cases considered here, 
the resulting values of the DR(n/p) are still far below the 
NSCL data. We notice that in the above calculation the 
isospin-dependent in-medium nucleon-nucleon scatte ring 
cross sections scaled by the nucleon effective mass [23j 
are used, and we have checked that the results do not 
change by much even if we use free nucleon-nucleon scat¬ 
tering cross sections. Thus, within the present IBUU11 
transport model, the variation of neither the symmetry 
energy nor the isospin effective mass splitting within their 
current uncertainty ranges is able to explain the exper¬ 
imental data. This situation clearly calls for possibly 
new mechanisms and explanations to resolve the neu¬ 
tron/proton ratio puzzle. It is thus interesting to note 
that several microscopic many-body theories and phe¬ 
nomenological models have predicted that the isospin 
dependence of short-range nucleon-nucleon correlations 
dominated by the tensor force can significantly reduce the 
kinetic part of the symmetry energy [5ll - l55| . The poten¬ 
tial part of the symmetry energy and thus the symmetry 
potential has to be enhanced accordingly to meet existing 
constraints on the symmetry energy at saturation den¬ 
sity. This effect has been shown to enhance significantly 
the double neutron/proton ratio of free nucleons in trans¬ 
port model calculations without considering the neutron- 
proton effective mass splitting nnm. A comprehensive 
study considering the isospin dependence of short-range 
nucleon-nucleon correlations and the neutron-proton ef¬ 
fective mass splitting in the IBUU11 model is planed. 

Trusting the data, the apparent success of the ImQMD 
model with the SLy4 interaction and the failure of the 
IBUU11 model with the ImMDI interaction in describ¬ 
ing the data requires further investigations. Generally 
speaking, different inputs and algorithms used in differ¬ 
ent transport models may lead to different predictions. 
We speculate, that the different handing of cluster for¬ 
mation might contribute appreciably to the difference be¬ 
tween the ImQMD and IBUU11 calculations. However, 
we honestly do not know at this time what are really 
the main causes because each model has several major 
inputs besides some technical differences. One possi¬ 
ble way out of this unfortunate situation is to conduct 
multi-observables versus multi-inputs covariance analy¬ 
ses, which have been successfully utilized in several other 
areas of nuclear physics recently (see Ref. [57j for a topic 


review). In the covariant analyses, the correlation matrix 
among all observables and input parameters as well as 
their uncertainties can be calculated consistently and si¬ 
multaneously. Such analyses with the IBUU11 model are 
underway. It will also be useful to perform such analyses 
with the ImQMD model. For example, the ImQMD cal¬ 
culations using the SLy4 and SkM* interactions predict 
an approximately 50% difference in the free n/p ratio at 
50 MeV/u. Where does this difference come from? It was 
attributed to the difference in the isospin effective mass 
splitting in Ref. [33; • In fact, at 5 = 0.2, m*/m n and 
m*/m p with SLy4 are respectively 0.68 and 0.71, while 
that with SkM* are respectively 0.82 and 0.76, with rel¬ 
ative effective mass splitting —3% for SLy4 and +6% for 
SkM*, smaller than that in the present calculation. How¬ 
ever, these are not the only differences between the two 
interactions. In particular, the isoscalar effective mass is 
different by about 14%, and the curvature of the sym¬ 
metry energy K sym differs by 30%. Without examin¬ 
ing sensitivities of observables to a particular input by 
varying it while fixing all others, it is hard to conclude 
what input is actually responsible for the observed change 
in any observable. Thus, multi-dimenstonal covariance 
analyses look promising although they are computation¬ 
ally extremely costly using transport models for nuclear 
reactions. 

In summary, within the IBUU11 transport model us¬ 
ing a newly improved isospin- and momentum-dependent 
interaction, we revisited but failed to resolve the neu¬ 
tron/proton ratio puzzle in heavy-ion collisions at inter¬ 
mediate energies. Nevertheless, some interesting physics 
and useful lessons are learned. We found that the 
neutron-proton effective mass splitting m* — m* indeed 
has a relatively stronger effect than the symmetry en¬ 
ergy E sym (p) and the assumption of to* < to* leads to 
a higher neutron/proton ratio of free nucleons and those 
in light clusters. Using the E sym (p) and to* — to* both 
within their current uncertainty ranges, with the IBUU11 
model and the ImMDI interaction we are unable to re¬ 
produce the recent NSCL/MSU double neutron/proton 
ratio data in central 124 Sn+ 124 Sn and 112 Sn+ 112 Sn col¬ 
lisions at 50 and 120 MeV/u. This situation clearly calls 
for new mechanisms to explain the puzzlingly high neu¬ 
tron/proton ratio observed in the experiments. Among 
the possible new physics origins, effects of the isospin- 
dependent short-range nucleon-nucleon correlation de¬ 
serve special attention. 
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